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ABSTRACT 

A  variety  of  experimental  tests  were  made  with  a  Gewa-T 
finned  surface  in  Freon-113  to  assist  in  understanding  the 
physical  mechanisms  which  occur  during  nucleate  pool  boiling 
from  this  type  of  surface. 

The  influence  of  the  channels  between  neighboring  T- 
shaped  fins  was  examined  with  the  use  of  four  aluminum  shrouds 
which  had  top  and  bottom  windows  of  various  apertures,  and 
which  were  fitted  tightly  over  the  Gewa-T  surface.   The  influ- 
ence of  the  fins  was  examined  by  progressively  machining  away 
the  fin  height  to  arrive  at  a  smooth  cylindrical  surface. 

The  application  of  the  shrouds  increased  the  heat  trans- 
fer coefficient  at  low  heat  fluxes  by  as  much  as  153  percent, 
but  decreased  the  heat  transfer  coefficient  at  high  heat 
fluxes  when  compared  to  the  unshrouded  surface.   The  addition 
of  the  "T-caps"  to  straight  fins  produced  the  most  significant 
improvement  in  heat  transfer  performance  when  compared  to  a 
smooth  tube.   Based  upon  these  results,  it  appears  that  the 
heat  transfer  performance  of  the  Gewa-T  surface  is  mainly 
influenced  by  the  effect  of  the  channels  created  by  neighboring 
T-shaped  fins. 
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I.   INTRODUCTION 

A.   BACKGROUND 

Nucleate  pool  boiling  is  characterized  by  bubble  forma- 
tion on  a  heated  surface  immersed  in  a  liquid  pool .   The 
bubble  nucleation  and  the  subsequent  bubble  separation  from 
the  surface  induce  fluid  motion,  or  mLxing,  near  the  surface 
which  enhances  heat  transfer  from  the  heated  surface  to  the 
liquid.   Techniques  to  enhance  boiling  heat  transfer  are 
being  vigorously  pursued  as  the  costs  of  labor,  material, 
and  energy  continue  to  escalate.   The  techniques  which  have 
received  the  greatest  attention  are  those  which  promote 
nucleation  as  a  result  of  surface  enhancement. 

The  purpose  of  surface  enhancement  is  to  provide  stable 
nucleation  sites  at  a  lower  wall  superheat  than  that  which 
is  required  for  boiling  from  the  natural  cavities  of  the 
surface,  and  thereby  produce  an  improvement  in  the  heat 
transfer  performance  over  that  of  the  unenhanced  surface. 
Improving  the  heat  transfer  performance  in  this  manner 
yields  a  corresponding  increase  in  the  heat  transfer  coeffi- 
cient due  to  the  lower  wall  superheat  required  for  boiling. 
Therefore,  the  enhanced  surface  requires  less  surface  area 
than  the  unenhanced  surface  for  the  same  amount  of  heat 
transfer.   Consequently,  a  savings  in  the  cost  of  material 
and  labor  is  achieved  by  reducing  the  surface  area  required. 


Enhanced  surfaces  have  been  classified  by  Yilmaz  and 

Westwater  [Ref.  1]  into  two  groups:   (1)  Porous  coatings,  and 

(2)  Integral  machined  surfaces.   Several  techniques  utilized 

in  the  production  of  porous  surfaces,  as  reported  by  Bergles 

and  Chyu  [Ref.  2]  are:   poor  welding,  sintering  or  brazing 

of  particles,  electrolytic  deposition,  flame  spraying,  bonding 

of  particles  by  plating,  galvanizing,  or  metallic  coating  of 

a  foam  substrate.   One  method  utilized  in  the  production  of 

integral  machined  surfaces  is  described  by  Webb  [Ref.  3] 

as  follows: 

State-of-the  art  concepts  employ  cold  metal  working  to 
form  a  high  area  density  of  reentrant  nucleation  sites, 
which  are  interconnected  below  the  surface.   Typically 
these  are  reentrant  grooves  or  tunnels,  rather  than 
discrete  cavities. 

Although  the  manner  in  which  these  two  classes  of  surfaces 
enhance  boiling  heat  transfer  is  not  fully  understood,  use- 
ful information  still  can  be  found  in  the  literature  about 
their  performance,  and  this  is  presented  below. 

The  heat  transfer  mechanisms  which  have  been  postulated 
for  boiling  from  porous  coated  surfaces  have  been  generalized 
by  Bergles  and  Chyu  [Ref.  4]  as  follows: 

...steady  vapor  formation  in  a  porous  matrix  is 
considered  to  be  primarily  internal.   There  are 
preferential  vapor  escape  channels,  the  liquid  is 
supplied  primarily  through  other  channels  surrounding 
the  vapor  channel.   Evaporation  occurs  at  the  mouths 
of  the  liquid  channels  and  on  adjacent  surfaces 
wetted  by  the  liquid.   A  low  AT  is  required  for 
evaporation  due  to  the  large  internal  surface  area. 
The  vapor  which  is  seen  in  the  pool  represents 
bubbling,  as  if  air  were  blown  into  the  pool  through 
a  perforated. plate,  rather  than  the  conventional 
nucleate  boiling  cycle  of  nucleation,  growth,  and 
departure . 
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In  the  case  of  the  machined  integral  surfaces,  similar 
mechanisms  may  occur . 

One  particular  integral  machined  surface  which  has  re- 
ceived attention  recently  is  the  Gewa-T  surface  manufactured 
by  Wieland-Werke,  AG.   During  19  81,  investigations  of  nucleate 
pool  boiling  from  a  Gewa-T  surface  were  conducted  by  Yilmaz 
and  Westwater  [Ref.  5],  Marto  and  Lepere  [Ref.  6],  and  Stephan 
and  Mitrovic  [Ref.  7] .   Yilmaz  and  Westwater  compared  the 
heat  transfer  performance  of  a  Gewa-T  surface  to  a  plain 

copper  tube  in  isopropyl  alcohol .   Their  results  show  that  at 

2 
a  heat  flux  of  4  0  kW/m  ,  the  heat  transfer  coefficient  of  the 

Gewa-T  surface  was  2.0  times  that  of  the  plain  tube.   Marto 

and  Lepere  compared  the  heat  transfer  performance  of  the  Gewa- 

T  surface  to  that  of  a  plain  copper  tube  in  Freon-113  and  FC- 

72.   Their  investigation  spanned  a  heat  flux  range  between 

2  2 

100  W/m  and  100  kW/m  .   Their  investigation  shows  that  for 

Freon-113  at  low  heat  fluxes  the  Gewa-T  surface  improves  the 

heat  transfer  coefficient  over  the  plain  surface  by  a  factor 

of  2.5,  and  at  high  heat  fluxes  by  a  factor  of  3.   Marto  and 

Lepere  also  determined  that  the  degree  of  superheat  required 

to  activate  nucleation  on  the  Gewa-T  surface  is  sensitive  to 

initial  surface  treatment  and  fluid  properties.   Stephan  and 

Mitrovic  compared  the  heat  transfer  performance  between  a 

Gewa-T  tube  and  a  standard  finned  Gewa  tube  in  R-ll,  R-12, 

R-22,  and  R-114 .   Their  investigation  covered  a  heat  flux 

2  2 

range  between  5  W/m  and  10  kW/m  .   The  results  of  their 


11 


investigation  show  the  heat  transfer  coefficients  of  the 
Gewa-T  surface  are  higher  than  those  of  a  tube  with  straight 
fins  of  the  same  diameter  and  same  fin  side  area.   Recently, 
Stephan  and  Mitrovic  [Ref .  8]  have  postulated  that  the  in- 
crease in  the  heat  transfer  coefficient  of  the  Gewa-T  tube 
over  the  standard  finned  Gewa  tube  occurs  in  the  following 
manner : 

A  single  vapor  bubble ,  formed  inside  the  T-shaped 
channel  of  the  Gewa-T  tube  after  break-off  flows 
around  the  tube  . . .  thus  keeping  contact  with  the 
tube  along  a  path  which  is  obviously  longer  than 
the  contact  length  around  the  Gewa  tube  with  plain 
fins. 

One  can  imagine  therefore,  that  the  bubbles  when 
moving  around  the  tube  sweep  off  other  forming 
bubbles  still  in  growth.   Thus  the  bubble  frequency 
and  also  the  heat  transfer  coefficients  increase. 
For  higher  heat  flux  densities  the  large  number  of 
bubbles  within  the  channel  grow  together  forming 
a  continuous  vapor  stream  in  the  core,  whereas  a 
thin  liquid  film  exists  near  the  wall. 

The  flow  pattern  around  the  tube  seems  to  be  of 
very  weak  influence  and  the  heat  transfer  process 
is  mainly  governed  by  the  bubble  formation  and  the 
two  phase  flow  inside  the  T-shaped  channel.   The 
top  part  of  the  fins,  which  were  rolled  to  form 
the  T,  practically  did  not  contribute  to  the  heat 
transfer  process,  because  even  at  maximum  heat 
fluxes  vapor  bubbles  were  not  produced  there. 

These  investigations  provide  an  indication  of  the  nucleate 

boiling  heat  transfer  performance  of  a  Gewa-T  surface,  using 

several  liquids  and  various  experimental  procedures.   Due  to 

the  sensitivity  to  initial  surface  treatment  and  fluid 

properties  exhibited  by  the  Gewa-T  surface,  it  is  not  possible 

at  present  to  generalize  the  heat  transfer  mechanisms  involved 

in  nucleate  pool  boiling  from  a  Gewa-T  surface  based  on  the 
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results  of  the  investigations  discussed.   However,  the  postu- 
late proposed  by  Stephan  and  Mitrovic  is  suited  to  experi- 
mental  testing. 

B.   THESIS  OBJECTIVES 

The  objectives  of  this  thesis  are  therefore  as  follows: 

1.  Design  and  construct  a  test  tube  from  a  Gewa-T  finned, 
solid,  copper  rod  which  will  permit  examination  of  nucleate 
pool  boiling  from  a  Gewa-T  surface  in  Freon-113. 

2 .  Gather  experimental  information  to  assist  in  under- 
standing the  physical  mechanisms  which  give  rise  to  increased 
heat  transfer  coefficients  with  this  type  of  surface. 


13 


II.   EXPERIMENTAL  DESIGN 

A.  REQUIREMENTS  GOVERNING  DESIGN 

The  purpose  of  the  experimental  apparatus  was  to  inves- 
tigate nucleate  pool  boiling  from  a  Gewa-T  surface  and 
attempt  to  ascertain  the  physical  mechanisms  which  give  rise 
to  this  phenomenon.   Consequently,  it  was  important  to  obtain 
reliable  measurements  of  those  parameters  which  may  provide 
an  indication  of  the  nucleate  boiling  mechanism.   These 
parameters  established  the  minimum  requirements  governing 
the  design  of  the  apparatus.   The  parameters  selected  for 
measurement  were  as  follows: 

1.  Ambient  pressure 

2 .  Surface  heat  flux 

3.  Test  section  wall  temperature  distribution 

4 .  Vapor  temperature 

5.  Fluid  bulk  temperature 

B.  TEST  APPARATUS 

The  test  apparatus  and  its  major  components  are  illus- 
trated in  Figure  1,  and  photographs  of  the  apparatus  are 
shown  in  Figures  2  and  3.   The  test  liquid  was  chosen  to  be 
Freon-113 . 

1 .   Boiler 

The  boiler  vessel  was  composed  of  a  thick-walled  pyrex 
glass  container  and  a  plexiglas  cover  fitted  with  a  rubber 
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o-ring.   The  vessel  had  an  outside  diameter  and  a  height  of 
30  5  mm,  and  these  dimensions  permitted  variations  in  test 
section  height  or  liquid  level .   The  vessel  was  made  of  glass 
to  provide  an  unobstructed  view  of  the  test  section  and  fluid 
during  the  collection  of  data,  and  it  was  situated  on  top  of 
an  electric  hot  plate  which  served  as  the  boiler  heating  ele- 
ment.  The  hot  plate  was  utilized  to  pre-boil  Freon-113  for 
degassing  prior  to  the  collection  of  data,  and  to  maintain 
the  fluid  at  saturation  temperature  during  operation.   Air  and 
noncondensable  gases  produced  from  the  degassing  process  were 
expelled  to  the  atmosphere  through  a  vent  identified  in 
Figure  1.   The  vent  to  the  atmosphere  permitted  the  boiler 
to  achieve  and  maintain  saturation  conditions  at  ambient 
pressure,  and  consequently  prevented  pressurization  of  the 
boiler . 

2 .   Condensers 

The  vapor  produced  from  boiling  was  condensed  with 
the  use  of  a  copper  primary  condenser  and  a  glass  secondary 
condenser.   The  condensers  were  water-cooled,  and  the  con- 
denser cooling  water  lines  were  connected  in  series.   The 
condensed  vapor  was  returned  to  the  liquid  bath  through  the 
action  of  gravity.   Consequently,  the  recovery  action  of 
the  condensers  minimized  the  escape  of  Freon-113  to  the 
atmosphere  as  a  result  of  venting  during  boiling,  and 
assisted  in  maintaining  the  conditions  of  saturation. 
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3.   Test  Section 

Figure  4  illustrates  a  cutaway  cross  sectional  view 
of  the  test  section.   The  outside  diameter  of  the  test  sec- 
tion was  21.2  mm  and  its  total  length  was  114.1  mm.   The  Gewa- 
T  surface  length  was  4  9.0  mm.   A  description  of  the  manufacture, 
assembly,  and  calibration  of  the  test  section  follows. 

A  spirally  finned,  solid,  copper  rod  designated  Gewa- 
T  13515.16,  manufactured  by  Wieland-Werke  AG,  was  utilized 
in  this  experiment.   Figure  5  illustrates  the  fin  profile  of 
this  surface.   The  characteristic  T-shaped  fins  were  formed 
by  rolling  a  standard  finned  Gewa  surface.   The  neighboring 
T-shaped  fins  of  this  surface  form  a  series  of  tightly  wound 
spiral  channels. 

a.   Machining  of  Test  Section 

A  radially-centered  hole  was  drilled  through  the 
length  of  the  solid  copper  rod  to  accommodate  a  cylindrical 
cartridge  heater.   Consequently,  a  tube  with  an  inside  diam- 
eter of  13.3  mm  was  manufactured  from  the  solid  rod.   Next, 
eight  longitudinal  holes  parallel  to  the  tube  axis  were 
drilled  in  the  tube  wall  up  to  the  axis  midpoint  at  a  radius 
of  8.0  mm.   All  eight  holes  had  a  constant  diameter  of  1.6 
mm.   Seven  of  the  eight  holes  were  drilled  at  intervals  of 
30  degrees  from  each  other  such  that  they  spanned  an  arc  of 
180  degrees.   The  eighth  hole  was  centered  at  an  interval  of 
90  degrees  from  either  of  the  two  nearest  holes.   Figure  6 
illustrates  the  angular  and  radial  position  of  the  holes  in 
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the  tube  wall.   The  Gewa-T  surface  was  then  machined  off  at 
each  end  of  the  tube,  producing  two  unenhanced  tube  ends  with 
an  outside  diameter  of  19.0  mm,  and  a  length  of  31.9  mm. 
Next,  the  inside  wall  of  the  unenhanced  ends  was  bored  out 
such  that  the  inside  diameter  of  the  ends  was  increased  to 
18.3  mm.   Consequently,  the  wall  thickness  of  the  unenhanced 
ends  was  reduced  to  0.35  mm,  creating  a  thin-walled  fin  at 
each  end.   This  action  was  taken  to  minimize  longitudinal 
conduction  losses  from  the  test  section. 

b.  Thermocouple  Assembly 

Eight  copper-constantan  thermocouples  were  manu- 
factured from  0.1  mm  diameter,  teflon- insulated  thermocouple 
wires.   The  average  thermocouple  bead  diameter  was  0.5  mm. 
These  thermocouples  were  inserted  into  eight  soft  copper 
capillary  tubes  whose  outside  and  inside  diameters  were  1.6  mm 
and  0.6  mm  respectively,  and  whose  outside  surfaces  were 
tinned  with  Eutectic  Corporation  157  solder.   The  thermocouples 
were  soldered  to  the  capillary  tube  end  through  which  they 
protruded  such  that  the  face  of  the  tube  end  served  as  a 
base  for  attachment. 

c.  Soldering  of  Thermocouples  and  Cartridge 
Heater  in  Test  Section 

A  cylindrical  cartridge  heater  was  radially  and 

axially  centered  in  the  test  section  with  an  aluminum  plug 

which  had  a  machined  cavity  on  its  interior  face  to  accept 

the  heater.   This  assembly  was  fitted  through  an  aluminum 

stand  and  placed  on  top  of  an  electric  hot  plate  which  had 
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been  heated  for  approximately  25  minutes.   The  hot  plate  was 
then  further  heated  to  a  temperature  of  approximately  350  °G 
in  order  to  bring  the  test  section  to  the  melting  temperature 
of  solder  which  was  218  °C.   Next,  a  bead  of  melted  solder 
was  dropped  into  each  of  the  eight  holes  bored  into  the  tube 
wall.   The  tinned  capillary  tubes  were  then  inserted  into 
these  eight  holes  and  soldered  into  place.   This  procedure  was 
immediately  followed  by  soldering  the  cartridge  heater  into 
place  using  the  same  Eutectic  Corporation  157  solder.   Upon 
completion  of  the  soldering  process,  the  assembly  was  removed 
from  the  hot  plate  and  quenched  under  flowing  tap  water. 

d.   Insulation  of  Test  Section  Ends 

After  the  test  section  was  cooled  to  ambient 
temperature  and  dried,  the  aluminum  plug  was  removed  and 
replaced  with  a  teflon  plug  of  the  same  dimensions.   The  end 
through  which  the  heater  and  thermocouple  wires  protruded 
was  insulated  with  poured  epoxy  resin.   The  thin-wall  surface 
at  the  teflon  plug  end  was  scribed  to  indicate  the  angular 
position  of  the  thermocouples. 

4 .   Calibration  of  Thermocouples 

Next,  the  0.1  mm  diameter  copper-constantan  thermo- 
couple wires  were  welded  to  0.3  mm  diameter  extension  wires, 
and  covered  with  heat  shrinkable  tubing.   The  thermocouple 
wires  were  wrapped  to  the  heater  wires  with  insulation  tape. 
The  test  assembly  was  placed  in  a  Rosemont  distilled  water 
constant  temperature  bath,  and  the  thermocouples  were  connected 
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to  a  Newport  digital  pyrometer  through  a  16-pole  thermocouple 
switch.   Two  separate  copper-constantan  thermocouples  were 
also  placed  in  the  bath  and  connected  to  the  same  Newport 
pyrometer.   The  test  section  thermocouples,  the  two  separate 
thermocouples,  and  the  Newport  digital  pyrometer  were  then 
calibrated  as  a  system  against  a  platinum  resistance  thermometer 
at  a  given  bath  temperature.   Each  thermocouple  temperature 
measurement  was  individually  read  with  the  use  of  the  thermo- 
couple switch  mentioned  above.   The  calibration  temperature 
range  was  between  39.23  °C  and  9  0.33  °C.   The  bath  temperature 
was  increased  in  increments  of  approximately  10  °C  and  then 
decreased  in  similar  increments,  except  for  an  initial  decrease 
of  6  °C. 

The  calibration  of  each  thermocouple  was  performed 
by  plotting  the  difference  in  temperature  between  the  thermo- 
couple and  the  platinum  resistance  thermometer  as  a  function 
of  thermocouple  temperature.   The  curve  which  best  fit  this 
data  was  a  straight  line  generated  from  a  least  squares  fit. 

Symbolically,  the  difference  in  temperature  between 
each  thermocouple  and  the  platinum  resistance  thermometer 
is: 


AT   =   TTC  -  Tp    (°C)  (1) 


where: 


'  c   =   thermocouple  temperature  (°C) , 


19 


Tp   =   platinum  resistance  thermometer 
temperature  (°C), 

AT   =  difference  in  temperature  between 

thermocouple  and  platinum  resistance 
thermometer  (°C). 


A  least  squares  fit  applied  to  the  data  yields: 


AT   =  A  TTC  +  B    (°C)  (2) 


where  A  and  B  are  constant  coefficients.   Therefore,  the 
calibrated  thermocouple  temperature  is  determined  from: 


TCAL   =   TTC  "  AT    (°C)  (3) 


where 


TCAL  =   calibrated  thermocouple  temperature  (°C) 


The  interval,  or  band,  centered  about  each  calibration  curve 
within  which  9  0%  of  the  data  were  located  determined  the 
uncertainty  in  temperature  measurement.   The  calibration 
curves  yielded  an  uncertainty  of  ±  0.15  °C. 

5 .   Installation  of  Test  Section  in  Apparatus 

At  the  completion  of  the  calibration  process,  the 
test  section  was  polished  and  cleaned  by  immersing  it  in  a 
10%  solution  of  nitric  acid  in  ethyl  alcohol  for  five  minutes 
After  properly  drying  the  test  section,  the  thermocouple 
and  heater  wires  located  at  one  end  were  inserted  into  a 
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copper  tube  fitted  with  a  Swagelok  fitting.   The  copper  tube 
had  a  gently  sloping  elbow  of  90  degrees  and  served  as  a 
support  arm  for  the  test  section.   The  Swagelok  fitting 
secured  the  test  section  to  the  support  arm. 

The  support  arm  and  the  two  calibrated  thermocouples 
were  secured  to  the  plexiglas  cover  on  the  boiler  vessel. 
One  of  the  two  thermocouples  was  immersed  in  the  freon  pool 
to  monitor  the  liquid  bulk  temperature,  and  the  other  was 
placed  above  the  boiling  fluid  to  measure  the  vapor  temperature 
6 .   Other  Instrumentation 

The  test  section  cartridge  heater  voltage  was  measured 
with  a  digital  voltmeter,  accurate  to  0.01  volts.   The  heater 
current  was  determined  by  measuring  the  voltage  drop  across 
a  2.031  ohm  precision  resistor  which  was  connected  in  series 
with  the  cartridge  heater. 

C .   SHROUDS 

Four  aluminum  shrouds  were  manufactured  with  an  outside 
diameter  of  22.2  mm  and  a  wall  thickness  of  0.4  5  mm.   Each 
shroud  measured  51.4  mm  in  length  and  had  two  diametrically 
opposed  windows  of  various  aperture  angles.   The  two  windows 
on  each  shroud  subtended  arcs  of  the  following  combinations: 
60  degrees  and  60  degrees,  30  degrees  and  30  degrees,  60 
degrees  and  30  degrees,  60  degrees  and  8.5  degrees.   The 
shrouds  were  slipped  over  the  test  section,  providing  a  tight 
mechanical  fit,  and  they  were  oriented  with  the  larger  aper- 
ture centered  vertically  at  the  top  of  the  test  section. 
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III.   EXPERIMENTAL  PROCEDURES 

The  work  of  Bergles  and  Chyu  [Ref .  9]  and  Marto  and 
Lepere  [Ref.  6]  show  the  dependence  of  boiling  incipience  on 
surface  past  history,  or  on  surface  treatment.   Therefore,  it 
was  necessary  to  establish  a  procedure  which  could  be  easily 
and  reliably  repeated,  and  which  would  permit  the  comparison 
of  data. 

A.   PROCEDURE  A 

Procedure  A  required  the  test  section  to  be  immersed  in 
the  liquid  pool  overnight.   The  liquid  level  was  maintained 
at  a  height  of  101.6  mm,  and  the  height  of  the  test  section 
from  the  base  of  the  glass  vessel  to  the  test  section  axis 
was  50.8  mm.   At  the  start  of  the  run,  the  hot  plate  voltage 
was  adjusted  to  the  maximum  setting  and  was  used  to  degas  the 
fluid  by  vigorously  boiling  the  pool  for  one  hour.   The  plate 
voltage  was  then  reduced  to  minimize  the  boiling  intensity 
and  maintain  the  pool  at  its  saturation  temperature.   Power 
was  then  supplied  to  the  test  section  heater  by  adjusting 
its  voltage.   The  heater  voltage  was  increased  in  increments 
between  2  and  30  volts,  and  then  decreased  in  increments  of 
approximately  10  volts.   The  system  was  allowed  to  stabilize 
for  5  minutes  at  each  power  setting  and  the  following  data 
were  recorded:   heater  voltage,  precision  resistor  voltage, 
vapor  temperature,  liquid  pool  temperature,  and  the  eight 
test  section  wall  temperatures. 
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B.  PROCEDURE  B 

Procedure  B  did  not  require  the  test  section  to  be  im- 
mersed in  the  pool  overnight.   Instead,  the  test  section  was 

2 
immersed  in  the  pool  and  subjected  to  pre-boiling  at  30  kW/m 

for  one  hour  using  the  cartridge  heater  while  the  hot  plate 
voltage  was  adjusted  to  the  maximum  setting.   Immediately 
after  this  one  hour  surface  aging  treatment,  the  power  to  the 
test  section  heater  and  the  plate  heater  was  secured,  and  both 
the  test  section  and  the  fluid  were  allowed  to  cool  for  one- 
half  hour.   The  plate  heater  voltage  was  then  adjusted  to 
the  saturation  temperature  setting,  and  the  cartridge  heater 
voltage  was  adjusted  in  various  increments  as  described  under 
Procedure  A.   At  each  power  setting,  the  same  data  as  de- 
scribe under  Procedure  A  were  recorded. 

C.  PROCEDURE  EMPLOYMENT  AND  TYPES  OF  RUNS 

Procedure  A  was  employed  only  once.   It  was  used  to  com- 
pare the  data  obtained  in  the  current  investigation  to  the 
data  of  Lepere  [Ref.  10]  and  to  observe  any  indication  of  the 
effect  test  section  diameter  may  have  on  the  heat  transfer 
coefficient.   In  this  investigation,  the  outside  diameter  of 
the  test  section  was  21.2  mm,  whereas  the  test  section  of 
Lepere  had  an  outside  diameter  of  17.9  mm. 

Procedure  B  was  utilized  in  all  the  other  data  runs. 
Two  data  runs  utilizing  Procedure  B  were  made  to  verify  the 
reproducibility  of  the  results  obtained.   For  this  reproduci- 
bility investigation,  the  test  section  height  was  50.8  mm 
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from  the  base  of  the  glass  vessel  to  the  test  section  axis, 
and  the  liquid  level  was  maintained  at  a  height  of  101.6  mm. 
In  two  other  reproducibility  tests,  the  test  section  was 
lowered  to  a  height  of  38.1  mm  in  one  data  run,  and  raised 
to  a  height  of  76.2  mm  in  the  other.   The  data  obtained  from 
these  two  runs  were  compared  to  the  data  in  which  the  test 
section  height  was  50.8  mm.   A  final  reproducibility  test 
was  performed  with  the  test  section  rotated  90  degrees,  and 
at  a  height  of  50.8  mm.   The  results  obtained  were  compared 
against  those  data  obtained  at  the  same  height  but  without 
the  rotation.   The  effects  on  nucleate  pool  boiling  due  to 
the  channels  formed  by  neighboring  T-shaped  fins,  and  the 
fin  geometry,  were  studied  utilizing  Procedure  B. 

The  channel  effect  was  investigated  in  four  separate 
data  runs  by  fitting  a  different  shroud  over  the  test  sec- 
tion for  each  run.   The  fin  effect  was  studied  by  progressively 
machining  away  the  fin  height  in  three  consecutive  data  runs. 
In  one  data  run  the  fin  "T-caps"  were  removed;  in  another 
data  run  the  fin  height  was  reduced  to  0.2  mm  above  the  base 
surface;  and  in  the  remaining  data  run,  the  fins  were  com- 
pletely removed,  and  the  tube  was  mechanically  polished  with 
emery   cloth  and  jeweler's  rouge  to  a  mirror  finish. 

D.   TEMPERATURE  CALCULATIONS 

The  saturation  temperature  was  determined  from  the  thermo- 
couple which  measured  the  vapor  temperature.   The  surface 
temperature  of  the  test  cylinder  wall  was  extrapolated  from 
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the  calibrated  temperature  measurements  of  the  eight  thermo- 
couples imbedded  in  the  tube  wall. 

The  surface  temperature  was  assumed  to  be  the  tempera- 
ture at  the  base  of  the  fins.   The  wall  temperature  was 
computed  by  taking  the  average  of  the  eight  calibrated  wall 
thermocouple  measurements,  and  then  a  correction  was  applied 
for  the  temperature  drop  across  the  tube  wall  due  to  conduc- 
tion in  the  test  cylinder.   Therefore,  the  surface  wall 
temperature  is: 

r2 

In  — 

Tb   "   TAVG  "  ^HT1    (°C)  (1) 


where: 


T,     =  wall  temperature  at  the  base  of  the 
fins  (°C), 

T   _   =   average  calibrated  thermocouple 
temperature  (°C), 

8 

TAVG   =   (  I      Tn)/8'°  (°C)  (2) 

n=l 


r,     =   radius  of  thermocouple  position  (m) , 

r2     =   radius  of  the  test  section  to  the  base 
of  the  fins  (m) , 

k     =   thermal  conductivity  or  copper  (W/m  -°C) , 

L     =   length  of  the  enhanced  surface  (m) , 

q     =   heat  transfer  rate  (W) . 
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The  heat  transfer  rate  was  determined  by  computing  the  power 
to  the  test  section  heater  in  the  following  manner: 

V     =   cartridge  heater  voltage  (volts) , 

vnt?c?   =   precision  resistor  voltage  drop  (volts)  , 

IvCtO 

R     =  resistance  of  precision  resistor  (ohms) , 


I   =   V^g/E   (amps)  (3) 


q   =   V  I    (W)  .  (4) 

E.   HEAT  FLUX  CALCULATIONS 

The  surface  heat  flux  was  determined  by  subtracting  the 
heat  losses  due  to  longitudinal  heat  conduction  through  each 
of  the  unenhanced  thin-walled  ends  from  the  power  supplied 
to  the  test  cylinder  heater,  and  then  dividing  the  differ- 
ence by  the  area  of  the  test  surface  at  the  base  of  the 
fins.   It  was  assumed  that  free  convection  from  a  horizontal 
cylinder  was  applicable  to  the  unenhanced  ends.   The  average 
heat  transfer  coefficient  of  the  unenhanced  ends  was  then 
determined  from  the  appropriate  Churchill  and  Chu  correla- 
tion [Ref.  11]  for  the  average  Nusselt  number: 

1/6 


o 


0  .387    RaD'  | 

|0*60    +    [1    i    (O.SSVPr)9/16]8'27  j 


NuD  =  [0.60     +    TTTTT—cry^n    ]  (5) 


10"5       <      RaQ         <      102. 
o 
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The  average  Nusselt  number  is: 


h  D 
NuD   =   -^  .  (6) 

o 


Therefore,  the  average  heat  transfer  coefficient  is: 

.      (  °'387    RaD/6  )2 

h      =      £-0.60+ °  (7) 

Do   |  [1    +    (0.559/Pr)9/16]8/27  ) 


where : 


k     =   thermal  conductivity  of  Freon-113 
(W/m  -  °C) , 

D_    =   outside  tube  diameter  at  unenhanced 
°       ends  (m), 

2 

v     =   kinematic  viscosity  (m  /s) , 

2 

a     =   thermal  diffusivity  (m  /s) 

Pr    =   -  (8) 

a 


ge(TW  -  TSAT)Do 
to     =   SAT  (9) 

O  V 


Tw    =   unenhanced  surface  wall  temperature  (°C), 
TSAT   =   fluid  saturation  temperature  (°C), 


Tw  -  TqA_   =   average  difference  between  surface 
wall  temperature  and  fluid 
saturation  temperature  (°C), 

2 
g     =   gravitational  acceleration  (m/s  ) , 

S     =   volumetric  thermal  expansion  coefficient 
(1/°C) 

=  T  It  (1/°c).  (10) 
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In  equation  (9) ,  however,  (T..  -  TqAT)  is  not  explicitly 
known . 

The  thin-walled  ends  were  modeled  as  rectangular  fins 
with  insulated  tips,  in  order  to  determine  the  heat  loss 
q£from  the  ends.   Therefore,  the  heat  loss  from  a  thin-walled 
fin  with  an  insulated  tip  at  each  end  of  the  test  section 
is  computed  from  the  following  equation  found  in  Incropera 
and  Dewitt  [Ref.  12]: 


=   [hPk  A]1/2  9,  tanh  mL     (W)         (11) 


where: 


P        =      tube   outside  wall   perimeter    (m) , 


P      =        ttDq         (m)  (12) 


k   =   thermal  conductivity  of  copper 
C      (W/m  -  °C) , 

2 
A   =   tube  cross-sectional  area  (m  ) 


A   =   J(D2  -  D2)    (m2)  (13) 

4   O     1 


D.   =   inside  tube  diameter  at  unenhanced  ends  (m) , 
t   =   fin  thickness  (m) , 


t   =   (D-D.)    (m)  (14) 

o     1 


L   =   Harper-Brown  length  correction  (m)  , 

L    =   L   +  %         (m)  (15) 

c       u    2 
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L   =   length  of  unenhanced  surface  (m) , 

m  =   (^-|)1/2    d/m)  (16) 

c 


It  was  assumed  that  the  temperature  at  the  base  of  the  thin- 
walled  fin  was  equal  to  the  temperature  at  the  base  of  the 
Gewa-T  fins.   Therefore: 


and 


9   =   TW  -  TSAT    (°C)  (17) 


9b   =   9 (0)   -   Tb  -  TSAT    (°C)  (18) 


cosh  m(L   -  x) 
__c 

cosh  mL 
c 


(19) 


The  average  difference  between  the  wall  temperature  and  the 
saturation  temperature  may  now  be  determined  from  the 
temperature  distribution  equation  as  follows: 


I.  cosh  m(L  -  x) 

-= c r— (°C)  (20) 

cosh  mL 


c 

,     c     cosh  m(L  -  x) 

9  ■  Z-l        %  coshL    dx    <°C>      (21» 

c  0  c 

eh 

tanh  mL     (°C) 


mL         c 


=   T   -  T       ( °C)  ( 22) 
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However,    m  is  a   function  of   R",    and   h   is   a   f unction  of   9  . 
This  yields   a   transcendental    equation   for  h: 


h    =    ~~  ;0.60  + 
o 


g66D^  1/6 

0.387( -) 

v  a 


) 


(W/m-°C)        (23) 


h-    k 
o 


0.60+0.387 


'geeb[tanh(^|)^2Lc]D^ 
c 

,h  P.  1/2  _ 

(FA}  Lc  V  a 

c 

[1  +   (0.559/Pr)9/16]8/27 


(W/nT  -  °C) 


(24) 


The  value  of  h  is  determined  by  an  iterative  technique. 
An  initial  guessed  value  for  h  is  substituted  into  the  right 
hand  side  of  equation  (24) ,  and  a  new  value  h1  is  computed. 
The  new  value  h*  is  now  subtracted  from  the  old  value.   If 
the  difference  is  less  than  .001,  then  the  value  of  h  has 
been  determined;  otherwise,  the  new  h1  is  substituted  into 
the  right  hand  side  of  the  above  equation  and  this  process 
continues  until  the  difference  is  less  than  .001. 

Once  h  is  computed,  the  heat  loss  may  be  determined  by 
substituting  the  value  of  h  into  equation  (11) .   After  the 
heat  loss  is  determined,  the  corrected  heat  transfer  rate 
from  the  test  section  is  computed  as  follows: 


^c   =   *  "  2<*£ 


(W) 


(25) 
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Therefore,  the  corrected  heat  flux  from  the  test  surface 
is: 


<*"   =  T^L  <W/m~»  (26) 
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IV.   RESULTS  AND  DISCUSSIONS 

Thirteen  experimental  runs  were  made.   The  data  obtained 
are  illustrated  in  the  figures  provided,  and  the  figures 
are  arranged  in  the  order  of  the  following  discussion. 

A.   COMPARISON  WITH  DATA  OF  LEPERE  [REF.  10] 

Procedure  A,  as  described  in  Chapter  III,  was  performed 
in  exactly  the  same  manner  as  reported  by  Lepere.   However, 
the  test  cylinder  used  in  this  investigation  had  eight  thermo- 
couples soldered  to  the  tube  wall   as  described  in  Chapter  II, 
whereas  the  test  cylinder  used  by  Lepere  had  four  equidistant 
thermocouples  soldered  to  the  surface  of  his  cylindrical 
cartridge  heater.   In  addition,  the  outside  diameter  of  the 
Gewa-T  surface  utilized  in  this  work  was  21.2  mm,  whereas 
the  outside  diameter  of  the  Gewa-T  tube  used  by  Lepere  was 
17.9  mm.   The  surface  geometry  of  the  two  test  cylinders  was 
identical  in  all  other  respects. 

Figure  7  illustrates  the  data  of  Lepere  plotted  together 
with  the  data  collected  in  this  investigation.   The  new  data 
conform  closely  to  the  Lepere  data,  and  a  hysteretic  behavior 
is  observed  in  both  sets  of  data.   At  low  heat  fluxes,  the 
agreement  is  excellent;  as  heat  flux  increases,  however,  the 
Lepere  data  exhibit  a  heat  transfer  coefficient  which  is  as 
much  as  1.4  times  the  heat  transfer  coefficient  of  the  Gewa- 
T  surface  utilized  in  this  investigation.   The  present  data 
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indicate  that  the  heat  transfer  coefficient  levels  off  in 
the  high  heat  flux  region  at  a  faster  rate  than  do  the  Lepere 
data.   While  a  quantitative  explanation  of  this  behavior  is 
not  provided,  a  qualitative  assessment  of  the  data  plotted 
supports  the  work  of  Cornwell,  Schuller,  and  Einarsson  [Ref . 
13] ,  which  examines  the  effect  of  tube  diameter  on  the  heat 
transfer  coefficient.   They  showed  that  increasing  the  outside 
tube  diameter  produces  a  decrease  in  the  heat  transfer 
coefficient. 

B.   REPRODUCIBILITY  TESTS 

The  data  obtained  from  the  reproducibility  tests  described 
in  Chapter  III  are  illustrated  in  Figures  8,  9,  10,  11,  and 
12.   Figure  8  is  selected  as  the  standard  plot  of  the  data 
collected  utilizing  procedure  B;  these  data  are  plotted  as 
a  solid  curve  in  all  subsequent  figures  for  the  purpose  of 
comparison.   Figure  9  compares  the  reproducibility  of  two 
sets  of  data  obtained  utilizing  Procedure  B  under  identical 
conditions.   The  agreement  is  excellent.   Figures  10  and  11 
compare  the  data  obtained  from  lowering  the  test  section 
height  and  raising  the  test  section  height  respectively, 
against  the  standard  data  depicted  by  a  solid  curve.   Once 
again,  the  agreement  is  excellent,  indicating  that  liquid 
level  in  the  pool  is  not  an  important  variable. 

The  final  test  of  reproducibility  is  plotted  in  Figure 
12.   In  this  case,  the  test  cylinder  was  rotated  9  0  degrees 
for  the  purpose  of  comparing  any  non-uniformities  in  the 
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peripheral  wall  temperature  distribution.   However,  it  was 
determined  that  the  uncertainty  in  temperature  measurement, 
which  was  ±  0.15  °C,  did  not  permit  a  good  comparison  since 
the  peripheral  temperature  variation  was  only  approximately 
±  0.30  °C.   Therefore,  the  temperature  data  reveal  the  neces- 
sity of  utilizing  a  more  precise  temperature  measuring  instru- 
ment in  order  to  make  a  comparison  possible.   Figure  12, 
however,  shows  that  the  data  obtained  from  rotating  the  test 
cylinder  9  0  degrees  are  in  excellent  agreement  with  the 
standard  data. 

It  is  apparent,  therefore,  that  data  may  be  reproduced 
utilizing  Procedure  B  with  excellent  agreement  to  the  stan- 
dard data.   There  is  no  recognizable  influence  due  to  the 
height  of  the  test  section  in  the  liquid,  nor  is  there  any 
apparent  influence  due  to  rotation  of  the  test  cylinder. 
The  figures  cited  support  the  confidence  placed  in  the  data 
obtained. 

C.   CHANNEL  EFFECT 

The  initial  study  of  the  channel  effect  on  nucleate  pool 
boiling  was  made  by  direct  physical  observation  of  the  Gewa- 
T  surface  during  operation.   At  low  heat  fluxes,  bubbles  were 
observed  to  flow  out  from  between  the  channels  located  at  the 
top  of  the  tube.   The  direction  of  the  bubble  flow  was  up, 
toward  the  surface  of  the  liquid.   As  the  heat  flux  was  in- 
creased, the  bubble  frequency  was  observed  to  increase,  and 
bubbles  of  about  1.5  times  the  diameter  of  those  leaving  from 
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the  top  of  the  tube  were  observed  flowing  out  from  channels 
on  the  lower  half  of  the  Gewa-T  surface.   In  addition,  bubbles 
of  nearly  3  times  the  size  of  those  at  the  top  were  observed 
forming  between  the  fin  tips  along  the  bottom  of  the  tube. 
These  large  bubbles  originated  from  smaller  bubbles  which 
were  pushed  out  through  the  gap  between  the  T-shaped  fins. 
The  bubbles  remained  situated  between  the  fins  and  continued 
to  grow  at  this  location  until  they  were  eventually  separated 
from  the  surface  by  the  action  of  local  convection  currents. 
At  high  heat  fluxes,  the  frequency  of  large  bubble  formation 
and  the  bubble  efflux  from  channels  along  the  bottom  of  the 
tube  were  so  great  that  these  bubbles  essentially  formed  a 
vapor  blanket  on  the  bottom  of  the  tube.   Similarly,  the  fre- 
quency of  formation  and  efflux  of  bubbles  from  the  remaining 
channels  on  the  lower  half  of  the  Gewa-T  tube,  and  from  the 
top  of  the  Gewa-T  tube  also  increased.   This  observation  has 
not  been  reported  in  the  literature  to  the  best  of  the  author's 
knowledge . 

The  observation  made  on  bubble  formation,  growth,  and 
efflux  from  the  channels  does  not  agree  with  the  description 
provided  by  Stephan  and  Mitrovic.   They  explain  that  bubble 
formation  and  growth  take  place  entirely  within  the  channels, 
and  bubbles  only  flow  out  from  channel  openings  located  toward 
the  top  of  the  Gewa-T  surface.   It  was  observed,  however,  that 
not  all  bubbles  flow  around  and  within  the  tube  channels. 
Bubbles  with  a  diameter  less  than  that  of  the  gap  between  the 
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fins  flow  through  the  gap  and  out  into  the  liquid  pool.   In 
addition,  at  high  heat  fluxes,  bubble  nucleation  was  observed 
to  occur  at  the  outer  surface  of  the  fins,  whereas  Stephan 
and  Mitrovic  reported  that  they  did  not  observe  this  condition 
during  their  investigation.   This  was  most  probably  because 
their  heat  fluxes  were  substantially  less  than  those  used 
during  the  current  investigation. 

As  heat  flux  is  increased,  bubble  generation  increases 
and  more  bubbles  exit  from  between  the  channels.   The  bubbles 
exiting  from  between  the  channels  at  the  bottom  of  the  tube 
disrupt  the  liquid  flow  upward  into  the  channels  and  eventually 
blanket  the  tube  surface  at  the  highest  heat  flux.   Conse- 
quently, the  heat  transfer  coefficient  levels  off.   Therefore, 
if  the  bubbles  can  be  constrained  to  flow  only  within  the 
channels,  an  increase  in  the  heat  transfer  coefficient  is  ex- 
pected.  The  shrouds  described  in  Chapter  III  were  utilized 
to  examine  this  effect. 

The  aluminum  body  of  the  shrouds  closed  the  channel  open- 
ings created  by  the  fin  gaps,  and  produced  channels  which  con- 
strained the  bubbles  and  liquid  to  flow  within  the  closed 
channels.   The  shroud  windows  permitted  liquid  to  flow  into 
the  channels  located  at  the  bottom  of  the  tube,  and  a  two- 
phase  mixture  of  liquid  and  vapor  to  flow  out  of  the  channels 
located  at  the  top  of  the  tube.   The  range  of  window  aperture 
combinations  manufactured  was  the  result  of  pursuing  a  com- 
bination which  would  maximize  the  heat  transfer  coefficient. 
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Figures  13,  14,  15,  and  16  illustrate  the  data  obtained 
utilizing  the  shrouds.   All  of  the  window  combinations  in- 
creased the  heat  transfer  coefficient  at  low  heat  fluxes. 
When  the  shrouds  were  fitted  over  the  Gewa-T  surface,  it  was 
observed  that  bubbles  which  form  and  grow  between  the  fin  tips 
at  the  bottom  of  the  tube  are  pushed  back  through  the  channel 
openings  created  by  the  gap  between  the  fins,  and  are  circu- 
lated within  the  closed  channels.   Consequently,  more  bubbles 
are  now  available  to  remove  other  bubbles  in  the  process  of 
growth  and  formation  as  these  bubbles  circulate  in  the  channel, 
giving  rise  to  an  improvement  in  the  heat  transfer  coefficient 
when  compared  to  an  unshrouded  surface.   However,  at  high 
heat  fluxes,  the  heat  transfer  coefficient  decreased.   This 
decrease  may  be  explained  by  the  high  frequency  of  bubbles 
exiting  from  the  channel  openings  at  the  bottom  of  the  tube. 
These  bubbles  disrupt  the  flow  of  liquid  into  the  channels, 
and  the  channel  surfaces  become  liquid  starved.   Consequently, 
the  thin  liquid  film  which  exists  at  high  heat  fluxes,  as 
described  by  Stephan  and  Mitrovic  [Ref .  7]  is  evaporated  and 
the  heat  transfer  coefficient  is  observed  to  decrease.   A 
comparison  of  the  improvement  in  the  heat  transfer  coefficient 

between  the  unshrouded  and  the  shrouded  Gewa-T  tube  follows. 

2 

At  a  heat  flux  of  10  kW/m  the  heat  transfer  coefficient 

was  increased  by  90  percent  with  the  60° -60°  shroud,  and  by 
111  percent  with  the  30° -30°  shroud  when  compared  to  the  tube 
without  the  shroud,  as  illustrated  in  Figures  13,  and  14, 
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respectively.   At  the  same  heat  flux,  the  heat  transfer 
coefficient  was  increased  by  123  percent  with  the  60°-30° 
shroud,  and  by  153  percent  with  the  6  0°-8.5°  shroud  when  com- 
pared to  the  unshrouded  tube,  as  illustrated  in  Figures  15, 

and  16,  respectively. 

2 
At  a  heat  flux  of  100  kW/m  ,  the  heat  transfer  coefficient 

was  observed  to  show  no  change  with  the  6  0° -6  0°  shroud,  but 
it  decreased  by  33  percent  with  the  30°- 30°  shroud  when  com- 
pared to  the  unshrouded  tube,  as  illustrated  in  Figures  13, 
and  14,  respectively.   At  the  same  heat  flux,  the  heat  trans- 
fer coefficient  decreased  by  20  percent  with  the  60°-30° 
shroud,  and  by  36  percent  with  the  60°-8.5°  shroud  when  com- 
pared to  the  unshrouded  tube  as  illustrated  in  Figures  15, 
and  16,  respectively.   Therefore,  at  low  heat  fluxes  the  bubble 
circulation  in  the  channels  is  enhanced  by  the  application 
of  the  shrouds  which  in  turn  enhances  heat  transfer.   The 
improved  circulation  is  due  to  the  buoyant  forces  which  pro- 
pel the  less  dense  vapor  bubbles  to  the  liquid  surface  at  the 
channel  opening.   However,  at  high  heat  fluxes,  the  window 
apertures  at  the  top  of  the  tube  were  not  sufficient  to  permit 
adequate  flow  of  vapor-liquid  out  of  the  channels.   Similarly, 
the  bubbles  exiting  from  the  channel  openings  at  the  bottom 
windows,  disrupted  the  liquid  flow  into  the  channels.   In 
addition,  the  shroud  body  itself  prevented  liquid  from  flowing 
into  those  channels  covered  by  the  shroud.   Consequently, 
the  heat  transfer  coefficient  was  observed  to  decrease  in 
the  high  heat  flux  region. 


38 


D.   FIN  EFFECT 

The  fin  effect  was  investigated  by  progressively  machin- 
ing away  the  fin  height.   Figures  17,  18,  and  19  illustrate 
the  data  obtained  after  removing  the  fin  "T-caps,"  the  fins, 
and  polishing  the  unfinned  tube  to  a  mirror  finish,  respec- 
tively.  A  comparison  of  the  heat  transfer  coefficients  ob- 
tained with  the  machined  geometries  was  made  at  a  low  heat 

2 
flux  of  10  kW/m  and  at  a  high  heat  flux  of  100  kW/m  .   The 

effect  of  surface  roughness  was  examined  by  comparing  the 
heat  transfer  performance  between  a  mirror  polished  smooth 
surface  with  that  of  an  unpolished  smooth  surface.   Surface 
roughness  produced  an  increase  in  the  heat  transfer  coeffi- 
cient of  90  percent  at  the  low  heat  flux  and  an  increase  of 
58  percent  at  the  high  heat  flux.   The  addition  of  the  fins 
to  a  smooth  tube  increased  the  heat  transfer  coefficient  by 
only  5  percent  at  the  low  heat  flux,  and  by  38  percent  at 
the  high  heat  flux.   This  implies  that  the  addition  of  fins 
did  not  greatly  enhance  the  heat  transfer.   The  addition  of 
the  "T-caps"  to  the  straight  fins  increased  the  heat  transfer 
coefficient  by  58  percent  at  the  low  heat  flux,  and  by  50 
percent  at  the  high  heat  flux.   Therefore,  the  heat  transfer 
performance  of  the  Gewa-T  surface  is  predominantly  influenced 
by  the  effect  of  the  channels  created  by  neighboring  T-shaped 
fins. 
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V.   CONCLUSIONS 

(1)  Bubbles  were  observed  to  burst  out  from  between  the 
fins  on  the  lower  half  of  the  tube,  and  nucleation  was  ob- 
served to  occur  on  the  outer  surface  of  the  fins.   Therefore, 
the  vapor-liquid  flow  path  is  not  confined  entirely  within 
the  channels  formed  by  the  T-shaped  fins  as  postulated  by 
Stephan  and  Mitrovic. 

(2)  The  effect  of  surface  roughness  was  examined  by 
comparing  the  data  of  the  mirror  polished  smooth  tube  to  that 
of  the  unpolished  smooth  tube.   Surface  roughness  increases 
the  heat  transfer  coefficient  by  90  percent  at  low  heat  fluxes 
and  by  58  percent  at  high  heat  fluxes. 

(3)  The  addition  of  fins  to  the  unpolished  smooth  tube 
increases  the  heat  transfer  coefficient  by  only  5  percent 
at  low  heat  fluxes  and  by  38  percent  at  high  heat  fluxes. 

(4)  The  addition  of  the  "T-caps"  to  the  straight  fins 
increases  the  heat  transfer  coefficient  by  58  percent  at  low 
heat  fluxes  and  by  50  percent  at  high  heat  fluxes. 

(5)  The  heat  transfer  coefficient  of  the  Gewa-T  surface 
may  be  increased  by  as  much  as  153  percent  by  using  shrouds. 

(6)  An  uncertainty  in  temperature  measure  of  ±  0.15  °C 
did  not  permit  comparison  of  the  tube  wall  temperature  pro- 
file.  A  temperature  measuring  instrument  with  a  precision 
of  at  least  0.01  °C  is  required  in  order  to  compare  wall 
temperature  profiles. 
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VI.   RECOMMENDATIONS 

The  following  recommendations  are  made  for  further 
study  of  the  nucleate  pool  boiling  mechanism: 

(1)  A  similar  experimental  study  should  be  made  with 
water  as  the  boiling  fluid.   A  comparison  can  then  be  made 
on  the  effect  of  bubble  diameter  between  fluids  with 
different  surface  tensions. 

(2)  A  study  should  be  made  on  the  influence  of  the  gap 
width  on  the  heat  transfer  coefficient  of  a  Gewa-T  tube. 

(3)  A  study  should  be  made  on  the  effect  of  surface 
irregularities  (e.g.,  indentations  along  the  fin  surface) 

on  external  nucleation  and  on  the  heat  transfer  coefficient. 

(4)  Testing  with  shrouds  should  be  continued  to  deter- 
mine the  optimum  shroud  which  maximizes  the  heat  transfer 
coefficient. 

(5)  A  study  should  be  made  to  compare  the  heat  transfer 
performance  between  a  standard-finned  Gewa  tube  and  a  Gewa-T 
tube  in  order  to  examine  the  effect  of  fin  geometry  on  tubes 
of  equal  surface  area,  and  identical  dimensions  in  all 
other  respects. 
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APPENDIX  A 
PROPERTIES  OF  FREON-113  [REF.  14] 


Boiling  Point 
(14.7  psia) 

Density,  116  (°F) 

Density,  118  (°F) 

Thermal  Conductivity 
(100  °F) 

Thermal  Conductivity 
(120  °F) 

Specific  Heat 
(100  °F) 

Specific  Heat 
(120  °F) 

Viscosity,  100  °F 

Viscosity,  120  °F 


117.63  (°F) ,  47.57  (°C) 

94.43  (lbm/ft3) 
94.26  (lbm/ft3) 
0.0419  (BTU/hr-ft-°F) 

0.0406  (BTU/hr-ft-°F) 

0.2  32  (BTU/lbm-°F) 

0.235  (BTU/lbm-°F) 

1.379  (lbm/ft-hr) 
1.229  (lbm/ft-hr) 
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APPENDIX  B 
SAMPLE  CALCULATIONS 

A.   WALL  TEMPERATURE  AT  THE  BASE  OF  THE  FINS 

The  surface  wall  temperature  at  the  base  of  the  fins 

is  computed  as  follows: 

r2 
In  — 

Tb  =  tavg  "  *{T71TY-)        (°c) 

c 
where  the  measured  parameters  are: 

V  =   20.74  (volts) 


VRES   =   1*36  (volts) 


R   =   2.0  31  (ohms) 


T±      =   53.4  (°C) 


T2   =   53.4  (°C) 


T3   =   53.3  (°C) 


T4   =   53.4  (°C) 


T5   =   53.3  (°C) 


Tc      =   53.3  (°C) 
6 


T?   =   53.3  (°C) 
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Tg   =   53.4  (°C) 


TSAT   =   47'6  (°C) 


D,   =   15.9  (mm) 


D2   =   19.3  (mm) 


L   =   49.0  (mm) 


k    =   383  (W/m-°c; 


Therefore, 


I   =  V^c/R      =   (1.36/2.031)   =   0.67  (amps! 


q   =   VI   =   (20.74) (0.67)   =   13.89  (W) 


8 

TAVG   =   (  I      V8*0)   =   53'4  (°C) 
n=l 

In  — 

rl 

Tb      TAVG  "  q  (  2  rr  L  k  } 

c 

"   53-4  "  13'89(2,(o"o4195),(383))    (°C) 


Tb   =   53.4  (°C)   (NOTE:   In  this  example,  the 

temperature  drop  is  negligible 

B.   HEAT  FLUX 

The  heat  flux  is  now  determined  at  the  given  heater 
voltage  setting  as  given  by: 
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"3"  =  TT^f-X   (w/n,2) 


where  the  parameters  required  are  obtained  as  follows 

k(100  °F)   =   0.0419  (BTU/hr-ft-°F) 
k(120  °F)   =   0.0406  (BTU/hr-ft-°F) 

Linear  interpolation  at  saturation  temperature  yields 


k(117.6  °F)   =   0.0419  +  (i^l-1-^0-)  (0  .  0406  -  0  .  0419) 

(BTU/hr-ft-°F) 


k(117.6  °F)   =   0.0408  (BTU/hr-ft-°F) ,  and 


converting  k  to  SI  units  yields: 


k(47.6  °C)   =   (0.0408) (1.7296)   =   0.0706  (W/ra-°C) 


Similarly: 


Cp(100  °F)   =   0.232  (BTU/lbm-°F) 


Cp(120  °F)   =   0.235  (BTU/lbm-°F) 


Cp(117.6  °F)   =   0.232  +  (H^^'Jjq0)  (0.235  -  0.232)  (BTU/lbm-°F) 
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Cp(117.6  °F)   =   0.235  (BTU/lbm-°F) , 


and  converting  C  to  SI  units  yields: 


Cp(47.6  °C)   =   (0.235) (4184)   =   983.2  ( joule/kg-°C) 


Similarly: 


y(100  °F)   =   1.379  (lbm/ft-hr) 


y(120  °F)   =   1.209  (lbm/ft-hr) 


y(117.6  °F)   =   1.379  +  (i^jfi-^gfi.)  (1.209  -  1.379)  (lbm/ft-hr) 


y(117.6  °F)   =   1.229    lbm-ft-hr, 


and  converting  y  to  SI  units: 


y(47.6  °C)   =   (1.229) (0.00041338)    (N-s/m2) 


y(47.6  °C)   =   5.080  x  lo"4    (N-s/m2). 


Similarly,  for  p: 


p(116    °F)       -      94.43    (lbm/ft3) 


p(118    °F)       =      94.26    (lbm/ft3) 
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p(117.6    °F)       =      94.43+    (i^i-Z^l)  (94  .  26  -  94  .  43)       lbm/ft3 


p(117.6    °F)       =      94.29    (lbm/ft3), 


and  converting  p  to  SI  units  yields: 


p(47.6  °C)   =   (94.29)  (  °-4536   )   kg/m3 

(0.3048)  J 


p(47.6  °C)   =   1510.41  (kg/m3). 


The  kinematic  viscosity  is  computed  as  follows: 

y              5.080     (10"9)  ,    oco     nrt-7,     ,    2,    . 

v      =      •=-     =      — 151Q    41 =  3.363     (10       )     (m  /s) 

The  thermal  diffusivity  is  computed  as  follows: 


k            0 .0706  .    _C/1  ,,  n-8>  ,  2  ,  . 

a      =   -p—  =   ~—  =   4.754  (10   )  (m  /s) 

P^P      (1510.41)  (9.817)  Q.0  ) 


The  volumetric  thermal  expansion  coefficient  is  computed  as 
follows: 


a  1  Ap_      ,     1 v  ,94.29  -  94.43, 

0        p  AT      l"  94.29M   117.6-116' 


6   =   0.00093   (A) 


1 

'F' 


converting  $  to  SI  units  yields: 
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B   =   (0.00093)(9/5)   =   0.0017   (^~) 


The  Prandtl  number  is  computed  as  follows 


Pr   =  1     =   3.363  (10"7) 
a      4.754  (10~8) 


Pr   =   7.0  7. 


The  remaining  parameters   are 


g      =      9.8       (m/s2) 


Dq      =       (.749)  (25.4)       =      19.0       (mm) 


P      =      ttDq      =      59.8       (mm) 


D.       =       (0.720)  (25.4)       =      18.3       (mm) 


A      =      tt(D2    -    D2)       =      8.202  x  10"5       (m2) 


t      =       (D      -   D.)       =      7.0  x  10    1      mm 


Lu      =       (1.255)  (25.4)       =      31.9       (mm) 


LC      =      Lu   +   I      =      31,9    +   ^~T~     =      32*3       (imn) 


)b      -      Tb   -   TSAT    CC)       =      53.4    -    47.6       (°C) 
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=   5.8   (°C) 


The  average  heat  transfer  coefficient  is  determined  from: 


h     =     £-  <  0.60  + 0.387- 


^8b[tanh(||1/2Lc3D3 

c 

,h  P,  1/2  r 
(kg>  Lc  V  * 


1/6 


[1  +   (0.559/Pr)9/16]8/27 


h      =      190.746       (W/m  -°C) 


The  heat  loss  may  now  be  computed  from 


[hPk  AJ1//2  9.  tanh  mL 


where : 


mL 


,h   P. 1/2 

lk  A; 
c 


=   0.60004 


(h  P  k  A)  1/2   =   0.61402 


Therefore,  the  heat  loss  from  each  unenhanced  end  is 


qz      =   2.137   (W) 


The  corrected  heat  transfer  rate  from  the  Gewa-T  surface  is 
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qc      =      q   -    2q        =      13.89    -    2(2.137)       =      9.62   W 


The  heat   flux  from  the  Gewa-T   surface    is: 


9    62  2 

q"      =      2_l£± _      =      3238       (W/iri   ) 

2.971    (10"l 
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APPENDIX  C 
UNCERTAINTY  ANALYSIS 

A.   UNCERTAINTY  IN  SURFACE  AREA 

In  order  to  determine  the  uncertainty  in  the  surface 
area,  it  is  necessary  to  first  determine  the  uncertainty  in 
the  diameter  at  the  base  of  the  fins,  and  the  uncertainty 
in  the  Gewa-T  surface  length. 


oi-   =  u)D   =   0.1   (mm) 


2 

A   =   it  D2  L   (m  ) 


where: 


D-   =   diameter  at  the  base  of  the  fins   (mm) 


L    =   length  of  the  Gewa-T  surface 


3A 

3*1  =  "  L 


3A        _ 

3L   "   *  D2 


Consequently,  the  uncertainty  in  the  surface  area  is 


uA   =   [(TTLajD  )2  +  (^D2ujl)2]1/2   (m2) 
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u)A      =       [(tt(49.0)  (0.1))2    +    (tt(19.3)  (0.1))2]1/2 


oj        =       [2  73.732]1//2      mm2 


a),      =      16.5      nun 

A 


The  uncertainty  in  the  surface  area  may  be  described  in  terms 
of  a  percentage  of  the  surface  area,  as  given  below: 


WA 

= 

WA 

A 

tt    D2    L 

WA 

16.5 

tt(19.3)  (49.0) 


WA 

-j-  =   0.0056,  or  0.56  percent. 


B.   UNCERTAINTY  IN  POWER 

The  uncertainty  in  power  is  determined  from  the  uncer- 
tainties in  the  test  section  cartridge  heater  voltage  and 
current.   The  measured  uncertainties  are  as  follows: 


ojt   =   0.01    (amps) 


oov  =   0.05    (volts) 


The  uncertainty  in  power  is  computed  as  follows 
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q   =   V  I    (W) 


j*  -   I    (amps) 


|^  =   V    (volts) 


2         2  1/2 
w    =   tdo)v)   +  (VWl)Z] 


2  2  1/2 

oo    =   [((0.67)  (0.05))^  +  (  (20.74)  (0.01)) *] 


oo    =   0.21 

q 


The  uncertainty  in  power  may  be  expressed  as  a  percentage  of 
the  power  as  given  below: 


CO  0J 

-%■     =   ,  ^  %<\      =   0.015,  or  1.5  percent 


C.   UNCERTAINTY  IN  AT 

The  uncertainty  in  AT  is  determined  from  the  uncertainty 
in  temperature  measurement.   The  uncertainty  in  temperature 
measurement  is: 


coT        =      oj  =      0.15         (°C) 

W  SAT 


AT       =       T       -    T 

W  SAT 
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|4S  =  1.0 

9Tw 


9AT-  =   -1.0 


9TSAT 


Therefore,  the  uncertainty  in  AT  is 


(joat   =   [  (ojt  )   +  (ajT    )  J 
W         SAT 


AT 


[(0.15)2  +  (0.15)2]1/2 


a)AT   =   0.2  °C 


Expressed  as  a  percentage  of  the  temperature  difference ,  the 
uncertainty  in  AT  is: 


WAT      0.2      n  n. 

-r-=r     =   -p-1-^  =   0.04,  or  4  percent. 


D.   UNCERTAINTY  IN  HEAT  FLUX 


q  "  2ql 

q   -  — s — 


5q"         <q-2q&) 
3A   *   -     A2 


3q"      1 

3q      A 
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3q"      _  2 

3q4      "      "A 

-q+2q«  2  12  2  2    1/2 


As   determined  earlier,    cu„    =   16.5   and  oj      =0.21.      In  addition, 

A  q 

it  was  assumed  that  cu    -  0 .  3q{  =  0.64.   Therefore,  the 


uncertainty  in  heat  flux  is: 


■•   „   [((   -9-62   fi)(16.5»2+  (    °-021  ,)2  +  (  2(0-64)  ?)2]1/2 
q         8.827(10"°)  2.971  (101      2.971  (10~J) 


w  "   =   4  37 

q 


The  uncertainty  in  heat  flux  as  a  percentage  of  the  total 
heat  flux  is: 


rr         4  17 

-gar     -   TTTft   =   0*135  or  13.5  percent 
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Figure   1.      Schematic   of  Test  Apparatus 
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Figure  3.   Photograph  of  Test  Section  in  Boiler 
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Figure  5.   Optical  Micrograph  of  Gewa-T  Fin  Profile, 


7X 
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Figure  6.   Thermocouple  Radial  and  Angular  Positions 
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Figure  7.   Procedure  A:   A  Comparison  of  Present  Data 
with  that  of  Lepere  [10] 
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Figure    8.      Procedure    B:      Standard   Data 
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Figure   9.       Reproducibility  Test:      Data    Under 
Identical   Conditions 
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Figure  10 .   Reproducibility  Test 
Lowered  to  3  8.1  mm 


Test  Section 
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Figure  11.   Reproducibility  Test:   Test  Section 
Raised  to  76.2  mm 
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Figure    12. 


Reproducibility  Test 
Rotated  90° 


Test  Section 
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Figure   13.      Application  of    60°-60°    Shroud 
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Figure   14.      Application  of    30°-30°    Shroud 
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Figure    15.      Application    of    6  0°- 30°    Shroud 
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Figure  16.   Application  of  60° -8. 5°  Shroud 


71 


6 
10     p 


i — i    i  nun 1 — i    mini 

:  0  INCRERSING 
•  t  DECREASING 
' —  DRTR  FROM  FIGURE    (8) 


T 1     I    M1IU 


10 


S  > 


13 


13 


13 


FREON-113 


J 1       1     I    I    I  111 1 I       I     '    I    "II  I  I! 


0.1 


1.0 


10.0 


100.0 


T   -   T         °C 
W         SRTy 


Figure  17.   "T-caps"  Removed,  Straight  Fins 
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Figure    18.      Fins    Removed,    Smooth  Tube 
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Figure   19 .      Mirror   Polished  Smooth  Tube 
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